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Abstract 

We quantify the important effect of strong final state interactions 
in the weak K — * 2tt amplitudes, using the measured tt-tt phase shifts 
with J = and / = 0,2. The main results of this analysis, with 
their implications for e'/e and the AI = 1/2 rule, have been already 
presented in a previous paper ||]. Here we provide a detailed formal 
derivation of those results and further discuss the Standard Model 
prediction of e'/e. 
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1 Introduction 



It is well known that, at centre-of-mass energies around the kaon mass, 
the strong S-wave tt-tt scattering generates a large phase-shift difference 
(<5g - 5g) (M 2 K ) = 45° ± 6° between the I = and I = 2 partial waves |]. 
In the usual description of K — > 27r decays, this effect is explicitly taken 
into account, through the following decomposition of the relevant isospin 
amplitudes with 1 = and 1 = 2: 

Aj = A[K -> (tttt)/] = A/ e i5 ° . (1.1) 



It is also known |3|-[12| that final state interactions (FSI) play an important 



role in the observed enhancement of the 1 = decay amplitude, A /A 2 ~ 
22.2. The presence of such a large phase-shift difference clearly signals a 
corresponding dispersive FSI effect in the moduli of the isospin amplitudes, 
because the real and imaginary parts are related by analyticity and unitarity. 
The size of the induced FSI correction can be roughly estimated from 



the available one-loop analyses of K — > 2n [ 12 - 14 1 in Chiral Perturbation 
Theory (xPT). At lowest order in the momentum expansion, 0(p 2 ) } the 
decay amplitudes do not contain any strong phase. Those phases originate in 
the final rescattering of the two pions and, therefore, are generated by chiral 
loops which are of higher order in momenta. Since the strong phases are quite 
large, one should expect large higher-order unitarity corrections. The one- 



loop calculations [O-Il] show in fact that the pion loop diagrams provide 
an important enhancement of the Aq amplitude, of about 40%. However, 
the phase-shift 5$ predicted by the one-loop calculation is still lower than 
its measured value, which indicates that a further enhancement could be 
expected at higher orders. 

Although the importance of FSI in K — > 2ir has been known for more than 
a decade, their impact on the direct CP- violating parameter e'/e has been 
overlooked in the so-called Standard Model predictions of this parameter, 
presented in refs. |15j and [fTB| . Not surprisingly, those predictions fail to 



reproduce the experimental measurements [[17 . 

The lattice investigations of kaon decay amplitudes have been only able, 
up to now, to compute the one-pion (k\TIas=i\K) matrix elements. In order 
to get the physical two-pion decay amplitudes, they rely on the lowest-order 
%PT relation between K — > ir and K — ► 27r, which, as mentioned before, 
does not include any FSI and underestimates the 1 = K — ► 2n amplitude 
by at least 40%. 



In refs. ||15|| , the large-iVc limit is used to fix the CP-violating K — > 2tx 



decay amplitudes. Since the strong phases 5j are zero at leading order in 
the 1/Nc expansion, the FSI enhancement has not been taken into account, 
either. Other approaches fl8|,|i"9| include some one- loop corrections and find 
larger values for the A amplitude. Although those are model-dependent 



1 



estimates, they provide an indication of the importance of higher-order pion- 
loop contributions. 

A proper way to account for the FSI effects has been addressed in ref. [[!]], 
where it has been shown that the strong rescattering of the two final pi- 
ons generates a large enhancement of e'/e. The resulting Standard Model 
prediction fllipOfl, 

Re(s'/e)\ SM = (17±6)xlQ-\ (1.2) 
is in good agreement with the present experimental world average f2l[ 



Re (e'/e) \ exp = (19.3 ± 2.4) x 10" 4 . (1.3) 



In the following we provide a detailed discussion of the approach advo- 
cated in ref. |1[ and further study its implications for e'/e and the AI = 1/2 
rule. The paper is organized as follows. In section |2| we formulate the Omnes 
problem for a general amplitude with two pions in the final state and derive 
its solution (for any number of subtractions). It is shown how the Omnes 
dispersive factor, solution of the Omnes problem, provides an all-order re- 
summation of the infrared chiral logarithms that contribute to FSI. This is 
a universal process-independent factor, which only depends on the quantum 
numbers (/ and J) of the final two-pion state. There is of course a poly- 
nomial (local) ambiguity, which encodes the process-dependent ultraviolet 
dynamics. 

To clarify the physics involved in the Omnes resummation, we present 
in section ^ explicit results for the scalar pion form factor. This quantity is 
known to two loops in the chiral expansion and, therefore, provides a sim- 
ple example where the power of our approach can be easily shown and the 
uncertainties quantified. The much more involved case of K — > 2tt transi- 
tions is discussed in section |], while section [5] presents the Standard Model 
prediction of e'/e. We conclude in section |6| with a few summarizing com- 
ments. We have relegated to the appendices the details on experimental tt-tt 
phase shifts, one-loop results for K — > 2ir and some remarks on recent 
literature on the subject. 



2 Omnes approach to FSI 

Let us consider a generic amplitude (or form factor) Aj(s), with two pions 
in the final state which have total angular momentum and isospin given by 
J and /, respectively, and invariant mass s = q 2 = (p\ + P2) 2 ■ As indicated 
by the name, FSI refer to the final dynamics of the two pions and not to the 
particular process leading to this final state. Therefore, we will not specify 
the physical amplitude and will look for a way to resum the strong n-n 
interactions to all orders in the chiral expansion. 
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Let us define the amplitude Aj(s) being analytic on the complex s plane 
except for a cut L = [4M%, oo) along the real positive s axis. For real 
values s < the amplitude is real; this implies that the values of the 

amplitude above and below the cut are complex conjugate of each other: 
A^ls + ie) = A T j(s - ie)*. Above the threshold, s > 4M^, A^s) has a 
discontinuity across the cut and develops an absorptive (imaginary) part. 

Cauchy's theorem implies that Aj(s) can be written as a dispersive inte- 
gral along the physical cut: 

. t i \ If , ImA^is) 

A\(s) = - dz J -—- + subtractions. (2.1) 

J it Jl z- s - ie 

The convergence of the dispersive integral is dictated by the specific form 
of the function Aj(s). Depending on the particular asymptotic behaviour 
of Aj(s) at the extremes of the cut L, a number of subtractions has to be 
performed to make the integral convergent. 

Let us further assume that Aj(s) corresponds to some weak or electro- 
magnetic transition, in the presence of strong interactions. Thus, above the 
cut, Aj(s+ie) = ((tt7t) j\0\i(q)) where O could be some effective (low-energy) 
electro weak Hamiltonian or a current. Working to first order in the small 
electroweak coupling, the unitarity condition allows then to write the imagi- 
nary part of Aj(s) as a sum over the contributions from all possible on-shell 
intermediate states which couple to the initial and the final state (properly 
normalized in momentum space): 

ImAfa + ie) = i£<Mj|rt|n) (n\0\i(q)) , (2.2) 

z n 

where T is the scattering T-operator. To derive this result, one makes use 
of Time-Reversal invariance; we will comment later on the proper way to 
bypass this assumption when analyzing CP-violating observables. 

Note that since ImAj is real, also the r.h.s. of eq. ( |2.2| ) is real. Below 



the first inelastic threshold, only the elastic channel contributes to the sum; 
one gets then: 

Im^ = (im A*j) = e' iS J sin^ Aj = e i& '> sin 8 T j A 1 / 

= sm5j \Aj\ = tan5jReAj, (2.3) 

implying that the phase of the decay amplitude Aj(s+ie) = \Aj(s) \ exp (^j) 
is equal to the phase of Tj, the im — > mi partial-wave scattering amplitude 
(Watson's theorem |22]|). Eq. ( |2.3| ) expresses the imaginary part of the am- 



plitude Aj(s) in terms of the amplitude itself, or its real part, and the titc 
phase-shift Sj(s). 
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Inserting eq. ( |2.3| ) in the dispersion relation ( [2.1Q , one obtains an integral 
equation for A(s) of the Omnes type, 



_ sq 7T J4M| [Z — Sq J -2 — S — IE 



(2.4) 



which has the well-known Omnes f|,^-|25| solution: 



A'M = Qj,n( s > s o) exp {l[ n (s, s )} , (2.5) 

where 

,j>,* )s <^r * _SM^ (2 . 6) 

7T J4M| (z — S ) Z — S — ie 



(n > 1) 



Qj,o( s > s o) = 1 and 

{n— 1 / „ \fc Jk 

S ife! ^ tog {^ (a) } 

(2.7) 

Strictly speaking, this equation is only valid below the first inelastic 
threshold (s < 16M^). However, the contributions from higher-mass in- 
termediate states are suppressed by phase space. The production of a larger 
number of meson pairs is also of higher order in the chiral expansion. 

We have written the most general result, for a given number of subtrac- 
tions n, performed at a generic subtraction point s outside the physical cut. 
The dispersive integral Ij n (s, s ) is uniquely determined up to a polynomial 
ambiguity (that does not produce any imaginary part of the amplitude), 
which depends on the number of subtractions and the subtraction point. 
This can be readily seen through the use of the following iterative formula 
for the real part of lj n (s, s ): 

Re 4,(5, s ) = Rel^s, s ) - (s - so)"" 1 lim Re I ^-i( s ^ ( {2 g) 

{ s — s o) 

where the second term on the r.h.s. depends on s only through the polynomial 
factor (s — so)™ -1 . The non-polynomial part of Ij n (s, Sq), containing the 
infrared chiral logarithms, does not have any dependence on the number of 
subtractions or the subtraction point. The polynomial ambiguity is of course 
canceled by the subtraction function Qj n (s, s ). 

Thus, the Omnes solution predicts the chiral logarithmic corrections in 
a universal way, independently of the number of subtractions or the sub- 
traction point, and provides their exponentiation to all orders in the chiral 
expansion. The polynomial ambiguity of Ij n {s, s ) and the subtraction func- 
tion Qj n (s,so) can be fixed, at a given order in the chiral expansion, by 
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matching the Omnes formula fl2.5|) with the %PT prediction of Aj(s). It 
remains a polynomial ambiguity at higher orders. 

A special case, which turns out to be relevant in the treatment of the 
weak K — > 2tt amplitudes, is the one where the amplitude Aj(s) has a zero 
of a given order p at some point s = (. In this case, once the zero is factorized 
through the relation Aj(s) — (s — Q p Aj(s), the dispersive representation of 
eq. (|2.5| ) is valid for the function Aj(s), so that 



A'As) = (s-CYA^s) = (s - C) p Qj, n ( s > s o) exp{/y S , So )}, (2.9) 



with Qj n (s, s ) the analogous of the expansion Qj n (s, s ) of eq. (|2^ ) for the 
function Aj(s) and Ij n (s, so) as defined in fl2.6|). 



3 The scalar pion form factor 

The scalar form factor of the pion is the simplest quantity where the Omnes 
problem can be solved |26],[27]] in order to resum final state interactions of a 
two-pion state with total angular momentum J = and isospin 1 = 0. It is 
defined by the matrix element of the SU(2) quark scalar density 



(n% P ')\uu + dd\7i k (p)) = 5 lk F*(t). 



(3.1) 



At low momentum transfer, %PT provides a systematic expansion of F$(t) 
in powers of t = {p ! — p) 2 and the light quark masses p8jp9(l : 



F* s (t) = F« s (0){l + g(t) + O(p 4 )} . 



(3.2) 



The value at t = coincides with the pion sigma term. It can be written as 
an expansion in powers of the light quark masses as follows: 



F£(0) 







+ 



d 



dm u dm<i 



Ml 



2B Q + 0(m 



(3.3) 



where Bq is a coupling of the lowest-order x?T Lagrangian which is related to 
the quark-antiquark vacuum condensate. The 0(p 2 ) correction g(t) contains 
contributions from one- loop diagrams and tree-level terms of the 0(p A ) %PT 
Lagrangian. It is given by [p8 , 29]: 



M, 



1 , 



Ml = 



9® = JUt) + l JKK(t) + T £J m (t)+4(Ll + 2Ll)(ri 



+ 



4(4tt) 



In 



Ml 



4(4tt) 



In 



M 



K 



Ml 



(3.4) 



with / w f w = 92.4 MeV. 
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Only two terms of the strong Lagrangian of 0(p 4 ) contribute to g(t). 
The /i dependence of their corresponding chiral couplings L r A (ii) and 
exactly cancels the one from the chiral loops appearing through the logarithm 
ln(// 2 /M^). At the standard reference value \i = M p , [L r 5 + 2L r ^{M p ) = 

(0.8 ±1.1) x lir 3 [gUHH. 

The functions J nn (i), Jkr^) and J w (t) are ultraviolet finite and, together 
with the logarithms, they are produced by the one-loop exchange of nir, KK 
and t]f] intermediate states. They have the form: 
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Below the first inelastic threshold, the absorptive part of the scalar form 
factor is generated by tttt exchange, through the one- loop function J^t). 
Thus, all non-analytic contributions originate from the final rescattering of 
the two pions and could be studied within the chiral SU{2) ® SU(2) frame- 



work [pqj , where the K and r\ modes are integrated out. In fact, for values 
of t such that t < 4M|>, 



J PP {t) 



It It 



2 



6 Ml 60 M 4 



(3.6) 



implying that, below the PP threshold, Jpp(t) has a very smooth behaviour 
and is strongly suppressed. In the case of the pion scalar form factor this 
means that at values of t <C 4M| r ,4M 2 the one-loop functions Jkr^) and 
J m {t) only give small analytic corrections, which are numerically negligible 
with respect to the local {L\ + 2L r A )(Mp) contribution. 

The two loop corrections to the pion form factor have been already com- 
puted |33]| . However, we prefer to keep the discussion at the one-loop level 



only, in order to make easier the comparison with K — > tttt where two-loop 
corrections are not yet available. 

Let us consider now the Omnes problem for Fg(t) and fix the subtraction 
polynomial performing a matching with the one- loop ;^PT result. We can 
write the Omnes solution in the form: 

F§(t) = n (t,t )-F£(t ) « Q o (t,t o )-F%(0) {l+g(t )} . (3.7) 

Thus, knowing the form factor at some low-energy subtraction point to, 
where the momentum expansion can be trusted, the Omnes factor Q (t,to) 
provides an evolution of the result to higher values of t, through the expo- 
nentiation of infrared effects related to FSI. The once-subtracted solution 
reads: 

fi?»(M„) = exp {till f _*L -JM-) s «<»((,(„) e *W . (3.8) 

I 7T JAM* Z — t Z — t — ie\ 
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We have split the integral into its real and imaginary part, making explicit 
that the phase of the Omnes factor is just the original phase-shift S^t). The 
Omnes exponential generates its corresponding dispersive factor 3£q (t, t ). 

The integral has been cut at the upper edge z, which represents the first 
inelastic threshold. Above z the representation ( |3.8|) is no longer valid and 
a coupled-channel analysis is required to solve the Omnes problem. From 
the behaviour of the S-wave 7T7t phase-shift Sq(z) (see appendix |A]), it is 
immediate to infer that the elastic integral evaluated up to z ~ 1 GeV 2 will 
slightly underestimate the exact result obtained with the inclusion of inelastic 
contributions. We shall discuss this point in more detail later. 

The solution for Fg(t) given in eq. (|3.7|) , being a physical quantity, 
must be independent of the subtraction point to, while the Omnes factor 
and the amplitude in front of it do depend on t . F° r illustrative pur- 
poses, we take t = M\ (the scale relevant for K — > tttt) and show in Ta- 
ble [I] the resulting value of \Fg (M^) / F$ (0)\ for different subtraction points 
t = 0, M 2 , 2M 2 , 3M 2 , AMI and m k- The u PP er limit of the integration 
range has been fixed at z = 1 GeV 2 . Our Omnes integral fl3.8| ) is not defined 
for t Q = M 2 K) because it lies above the threshold of the non-analyticity cut; 
however, by its definition in eq. ( |3.7|) , Q(Mx,to = M^) = 1 holds. 



to 


g(t ) 


3C(M!,t Q ) 


^\Ml,to) 


l^(^)/^(o)l 


n = 1 


n = 2 








(1.23) 1.45 


(1.55) 1.56 


(1.23) 1.45 


(1.55) 1.56 


Ml 


0.042 


(1.21) 1.40 


(1.47) 1.44 


(1.26) 1.46 


(1.53) 1.50 


2Ml 


0.091 


(1.17) 1.34 


(1.38) 1.31 


(1.28) 1.46 


(1.51) 1.43 


3MI 


0.15 


(1.12) 1.26 


(1.26) 1.13 


(1.29) 1.45 


(1.45) 1.30 


AMI 


0.26 


(1.03) 1.11 




(1.30) 1.40 




Ml 


0.54- 0.46 i 


= 1 


= 1 


1.61 


1.61 



Table 1: The one-loop function g(t ), the Omnes factor W ,(t, to) and the 
modulus of Fg(t)/Fg(to) are shown at t — M\ for different values of the 
subtraction point to € [0, Mjy and for n = 1, 2 subtractions. At a given to, 
the first value (within brackets) is obtained with the 0(p 2 ) %PT prediction 
for <5q, while the fit to the experimental phase-shift data has been used in the 
second one. The integrals have been cut at z = 1 GeV 2 . 

The dominant contributions to g(to) come from the logarithms and the 
7T7T one-loop function J n7T (t). The corrections from J^k(^o) and J m (to) stay 
within 1% of the total one-loop correction at all non-zero subtraction points, 
while the local [L r 5 + 2L r 4 ](M p ) term gives a contribution smaller than 10%. 
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The xPT calculation of Fg(to) is obviously better at lower values of to where 
the one-loop correction g(t ) is smaller. At to = 4M 2 a sizeable 26% effect is 
already found, while at t = M\ the correction is so large than one should 
worry about higher-order contributions. The Omnes exponential allows us 
to predict Fg(M^) in a much more reliable way, through the evolution of 
safer results at lower t values. 

The resulting values for the scalar form factor remain very stable within 
the whole range of subtraction points. Increasing the value of t one just 
moves higher-order x?T corrections from the Omnes factor to the amplitude 
Fg(to). At to = Ml the one-loop corrections are still almost zero, while the 
Omnes factor contains all the higher-order effects. At the highest possible 
subtraction point to = 4M 2 , i.e. the threshold of the non-analyticity cut, 
the bulk of the higher-order corrections has been moved to Fg(t ), while 
the Omnes factor approaches one, the value that it assumes at t = M|-, by 
construction. 

Together with the more accurate results obtained with the experimental 
phase-shifts, we have shown in Table [I] (within brackets) the corresponding 
numerical values using the lowest-order x?T prediction for 5q(z) in eq. (|A.4j) . 
As expected, the 0(p 2 ) x?T approximation to Sq(z) underestimates the dis- 
persive integral and, therefore, the true results, which are obtained with the 
experimental phase-shifts. 

As shown in Table |2], the dispersive correction factors ^^(M^, i ) in- 
crease with increasing values of the upper integration limit z. This clearly 
indicates that we are underestimating the FSI effect. However, one cannot 
trust the numerical results obtained for z > 1 GeV 2 because a coupled- 
channel analysis is required above the inelastic threshold. 





^\M 2 Kl t ) 


^o 2) (M|,t ) 


z= 1 


z = 2 


z = 3 


z = 1 


z = 2 


z = 3 


t = 


1.45 


1.58 


1.62 


1.56 


1.58 


1.59 


to = Ml 


1.40 


1.51 


1.55 


1.44 


1.46 


1.47 


to = 2Ml 


1.34 


1.44 


1.47 


1.31 


1.32 


1.33 


to = 3M 2 


1.26 


1.35 


1.38 


1.13 


1.14 


1.14 


to = 4M 2 


1.11 


1.19 


1.21 









Table 2: Dependence of 9f$ n) {M 2 K , t ) on the upper edge z (in GeV 2 units) of 
the dispersive integral, for various choices of t and n = 1, 2. The fit to the 
experimental data for 5$ has been used. 

We can suppress the sensitivity to the high integration range by using 
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a twice-subtracted dispersion relation. The corresponding Omnes factor is 
given by: 



«?>(*,*,) = exp( (i - (o)T ^) + (*-<») 2 r «w 



1 + 0(*o) J4Af| (2 - t ) 2 z-t-iej 

= M?J(t,t )e«8(*). (3.9) 

Table |2| shows that with n = 2 the numerical results remain indeed stable 
under variations of z. Moreover, as shown in Table [l], the 0(p 2 ) approxima- 
tion to 5q(z) works now much better, giving results in good agreement with 
the ones obtained from the experimental phase-shifts. 

Notice, that Q^\t,t ) is not defined at t Q = 4M^, because the derivative 
g'(to) has a discontinuity at the threshold of the physical cut. Since F$(t) 
is an analytic function in the cut s-plane, its Taylor expansion around the 
subtraction point to (< 4M^) has a convergence radius |to — 4M^|, which 
becomes zero at to = 4M%. Thus, subtraction points close to this threshold 
singularity [0 should be avoided |35| . 



Since the derivative g'(to) has been fixed at the one-loop level only (i.e. 
has been estimated at the lowest non-trivial order), there is a correspond- 
ing uncertainty which gets somehow increased by its exponentiation. This 
explains why the predicted values of \Fg (M^) / Fg (Q)\ in Table [I] are less 
stable for n = 2 than for n = 1 under changes of the subtraction point 
to- A twice-subtracted Omnes solution requires a more precise knowledge 
of the subtraction function. With n = 2, the scalar form factor slightly de- 
creases for increasing values of t , in the same way as 3t^\t, t ), the tree-level 
once-subtracted solution, does. This could be easily improved by using the 
available two- loop xPT results |$3| . Nevertheless, since chiral corrections are 



smaller at lower values of to, we can safely conclude that the true value of 
\Fg (M^/Fg (0)| is between 1.5 and 1.6. Taking the experimental phase-shift 
uncertainties into account, we finally get: 

\ f s( M k)/ F s{°)\ = 1-55 ±0.10. (3.10) 

The naive one- loop xPT prediction, \F^(M^)/F^(0)\ = 1.61, turns out to 
be within the la range of our final result ( |3.10| ). The advantage of using the 
dispersive Omnes resummation is that one can pin down the true value with 
an acceptable accuracy (7%), in spite of having a 60% one- loop correction. 



4 K —> 7T7T amplitudes 

The analysis of K — ► tttt transition amplitudes is technically more compli- 
cated. Since the electroweak scale M\y, where the short-distance quark tran- 
sition takes place, is much larger than the long-distance hadronic scale, there 
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are large short-distance logarithmic contributions which can be summed up 
using the Operator Product Expansion [36] and the renormalization group. 
One gets an effective AS = 1 Lagrangian, defined in the three-flavour the- 
ory |37|, |3~8|1 , 

£ e A / =1 = ~^V ud v: s Qi(y) , (4.1) 



which is a sum of local four-fermion operators Qi, constructed with the light 
degrees of freedom, modulated by Wilson coefficients Ciiv) which are func- 
tions of the heavy masses M w , M z , m t , m h and m c that have been integrated 
out. The overall renormalization scale v separates the short- [M > v) and 
long- (m < v) distance contributions, which are contained in Ciiv) and 
Qi, respectively. The physical amplitudes are of course independent of v; 
thus, the explicit scale (and scheme) dependence of the Wilson coefficients 
should cancel exactly with the corresponding dependence of the Qi matrix 
elements between on-shell states. We have explicitly factored out the Fermi 
coupling Gp and the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements 
Vij containing the usual Cabibbo suppression of K decays. 

Our knowledge of AS = 1 transitions has improved qualitatively in recent 
years, thanks to the completion of the next-to-leading logarithmic order cal- 



culation of the Wilson coefficients |39],^0[]. All gluonic corrections of 0(a"t n ) 
and 0(a™ +1 t n ) are already known, where t = In (Mi/M 2 ) refers to the log- 
arithm of any ratio of heavy mass scales M±,M 2 > v. Moreover, the full 
m t /My/ dependence (at lowest order in a s ) has been taken into account. In 
order to predict physical amplitudes, however, one is still confronted with the 
calculation of hadronic matrix elements of the four-quark operators. This is 
a very difficult problem, which so far remains unsolved. 

The chiral symmetry properties of the effective Lagrangian ( |4.1| ) deter- 
mine its corresponding xPT realization, in terms of the QCD Goldstone 
bosons 



V 



K- 



7/ 



7T 



2" 



hi 



K+ 
K° 



K c 



(4.2) 



parametrized through the exponential U = exp(^/2i& / f) . At a given order 
in the momentum expansion, chiral symmetry fixes the allowed chiral opera- 
tors and, therefore, the structure of the physical weak amplitudes. The only 
remaining problem is the calculation of the chiral couplings from the effective 
short-distance Lagrangian. 

The effect of strangeness-changing non-leptonic weak interactions with 
AS = 1 is incorporated |41| in the low-energy chiral theory as a perturba- 
tion to the strong effective Lagrangian. At lowest order, the most general 
effective bosonic Lagrangian, with the same SU (3) SU (3) r transformation 
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properties as the short-distance Lagrangian (O 



), contains three terms: 



r AS=l 
*~-2 



'13 



-^V ud V: s S^g 8 f (\L,L») + g 27 f + h, 21 L{ 

+ e 2 fg EW (\U^QU)). (4.3) 



The flavour-matrix operator = —iWD^U represents the octet of V — A 
currents at lowest order in derivatives, Q = diag(|, — |, — |) is the quark 
charge matrix, A = (A 6 — i\ 7 )/2 projects onto the s — > d transition [A^- = 
^3^2] and (A) denotes the flavour trace of A. 

The chiral couplings g 8 and g 27 measure the strength of the two parts 
of the effective Lagrangian ( |4.1| ) transforming as (8l,1r) an d (27^,1^), re- 
spectively, under chiral rotations. Chiral symmetry forces the lowest-order 
Lagrangian to contain at least two derivatives (Goldstone bosons are free 
particles at zero momenta). In the presence of electroweak interactions, 
however, the explicit breaking of chiral symmetry generated by the quark 
charge matrix Q induces the 0(p°) operator (XWQU) J42],[43|], transforming 
as (8l,&r) under the chiral group. In the usual chiral counting e 2 ~ 0{p 2 ) 
and, therefore, the g E w term is also of order p 2 . One additional term E31 
proportional to the quark mass matrix, which transforms as (8^,1^), has 
not been written since it does not contribute to the physical K — > tttt ma- 
trix elements [12|, |I~3|. |45[ 



At next-to-leading order in the chiral expansion, 



i.e. 0(p ), a set of additional weak counterterms will contribute fl2|, |13|. [46 



together with the strong chiral operators Li introduced in 29 



The Lagrangian (|4.3| ) gives the lowest-order contribution to the K — > 2ir 
matrix elements. At generic values of the squared centre-of-mass energy 
s — (Pvri + Pw2) 2 , the 7 = 0,2 amplitudes are given by 



72 
G 



V ud V: s V2f {{98 + I 927) (s - Ml) - \fe 2 9Ew\ 
2 



A (s) = 

A 2 (s) = -^=V ud V: s ^f {5g 27 (s-Mt) - 3 / 2 e W} . 
We have made the usual isospin decomposition: 



(4.4) 



A + -^=A 2 . 



A[K° -> 7r + 7r] 

r°7T ] = An - V2A 

-0l 



A[K° 
A[K + - 



(4.5) 



7T 7T 



3 , 

- 2 A 2 , 



where the amplitudes Ai = Aj exp {i^} contain the strong phase-shifts, 
which are zero at tree level. 
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For the discussion of the CP-conserving amplitudes we will neglect^ the 
tiny electroweak correction proportional to e 2 gEw- Taking the measured 
phase-shifts into account, eqs. ( |4.4[ ) allow us to extract the lowest-order weak 
couplings from the experimental information on K — > 2tt decays |)0) : 



+ l 927 



5.1, \g„/ga\~ 1/18. (4.6) 



The huge difference between these two couplings shows the well-known en- 
hancement of the octet |AJ| = 1/2 transitions. 

Let us now apply the Omnes procedure to the AS = 1 decay amplitudes. 
This is more subtle than for the scalar pion form factor, because we need 
to consider an off-shell kaon of mass squared s = (p^i + P^Y instead of a 
physical momentum transfer s. Since we are just studying the corrections 
induced by FSI between the two pions, the kaon can be formally considered 
as an external source, provided all SU(3) symmetry constraints are satisfied. 
As we saw explicitly for the scalar form factor, the FSI corrections that are 
summed up through the Omnes exponential are actually an SU(2) effect, 
generated by pion loops. Intuitively, we are just correcting a local weak 
K — > 7T7T transition with a chain of pion-loop bubbles, incorporating the 
strong 7T7T — ► 7T7T rescattering to all orders in the chiral expansion. 

In the absence of e 2 gEw corrections, the tree-level isospin amplitudes 
have a zero at s — M%, because the on-shell amplitudes should vanish in 
the SU(3) limit [4~3, 4~7,ET, 52]. We must take this important constraint into 



account, when making the Omnes summation of FSI effects, factorizing the 
zero explicitly as indicated in eq. ( |2.9| ). This is what was done in ref. [|I|], 
using a once-subtracted Omnes factor, to evolve the tree-level x?T results 
from s = Ml to the physical point s = M\. 

At higher orders in xPT there are small corrections proportional to [M\ — 
M%) instead of (s — M%), which originate in the explicit breaking of chiral 
symmetry provided by the quark mass matrix. According to the general one- 
loop analysis presented in appendix ||, those tiny effects can be neglected to 
a very good approximation. However, there is no need to do it. In full 
generality, the isospin amplitudes can be decomposed as 

Ms) = a T (s) (s - Ml) + 5~ ai {s) (Ml - Ml) , (4.7) 

where Sai(s) is zero at lowest ordeiQ. Since there is a single strong phase, 
for a given isospin, the unitarity relation ( p.3|) is also valid for the individual 
functions ai(s) and 6ai(s). Therefore, the Omnes problem can be solved 

3 A general analysis of isospin breaking and electromagnetic corrections to K — > 2ir 
transitions is presented elsewhere 



4 To make the decomposition (4.7) unique, we require the function Sai(s) to depend 
on s only logarithmically. 
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separately for the two pieces. Combining them, we can write our result for 
the physical on-shell amplitude in the simpler form: 

Ai = MM 2 K ) = (M 2 K -Ml) ai (M 2 K ) 

= (M| - Ml) Q!(M 2 K , s ) aj(so) (4.8) 

= {Ml - M 2 ) MM 2 K , s ) ai (s ) , 

where a/(s) = a/(s) + 5di(s). 

The once-subtracted Omnes factor fij (Mj^, so) is universal, i.e. the 
same for f(s) = Fg(s), or 5ai(s), because it only depends on the 

strong phase-shift ^o(s) [see eq. ( 3.81) 1 . This is no longer true with two sub- 
tractions [eq. Q3.9D 1) because fi^^So) contains an explicit dependence on 
f'{so)/f(so)- Nevertheless, given the smallness of the non-leading Sai(s) 
contribution, it is a very good numerical approximation to take also a global 
Omnes exponential for Of(s) in the twice-subtracted case. 

Let us define ao(s) = a 8 \s) + a 27 \s), thus separating the (8l, 1r) and 
(27 l, 1r) contributions to the isoscalar amplitude. The complete one-loop 
%PT results for the different decay amplitudes are given in appendix |B[ 
Their s dependences can be written in a rather transparent way: 

„(8)^ _ 



= 4 8) (o) {i + gl?\s)+0(p 4 )} , 

= 4 27) (0) {l+9o [27 \s)+0(p 4 )} , (4.9) 
a 2 (s) = a 2 (0) {l+g 2 (s)+0(p 4 )} , 



,(27) 



where 



(4.10) 



9T\s) 




(4.11) 



(4.12) 



The corrections coming from the Sai terms, included in these results, 
are very small. Denoting by gi(s) the corresponding functions for the un- 
corrected aj amplitudes, the differences Ag T (s) = gi(s) — gi(s) only get 
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contributions from the KK and 7777 loop functions which, as shown in sec- 
tion |, are numerically suppressed at low values of s. Moreover, they get 
an additional suppression factor M 2 (see appendix B). Since the KK and 
7777 intermediate states cannot give rise to I = 2, one gets Ag 2 (s) = (5a 2 
does not have any s dependence at this order), while the isoscalar differences 
are given by Ag ( *\s) = M^[9J KK (s) + 8J m (s)]/(3Qf 2 ) and Ag ( * 7 \s) = 
Ml [3J KK (s) - 4J w (s)]/(2/ 2 ). Even at s = M%, these differences are com- 
pletely negligible: Agf\M] c ) ~ 1 x KT 3 and Ag^ 27) (M%) 2 x 10" 4 . 

Notice the strong similarity with the scalar form factor result in eq. ( |3.4j) . 
The isoscalar functions g$ (s) and g^ 7 \s) get exactly the same J^r(s) con- 
tribution than the scalar form factor function g(s), while the corresponding 
contribution to ^( s ) has opposite sign. The polynomial factors in front of the 
7T7T loop function, (s — M 2 /2) and (M 2 — s/2) for 1 = and 2 respectively, 
clearly identify the corresponding lowest-order xPT phase-shifts, given in 
eq. fliU) [Ini J n Js) = 9(s - 4M 2 ) <tM/(1&k)]. 

In addition, at /i 2 = M|, one recognizes the same local infrared log- 
arithmic enhancement of all isoscalar g functions, In (M| r /M 2 )/(47r) 2 ; the 
corresponding factor in g 2 decreases the 1 = 2 amplitude. For arbitrary val- 
ues of the chiral scale /x, this logarithmic correction is split in In (/x 2 /M 2 ) and 
In (M|-/M 2 ) terms, which are slightly different for the different g functions. 
However, all isoscalar g functions contain exactly the same infrared InM 2 
contribution. 

Thus, the s dependence of the weak decay amplitudes is indeed dominated 
by infrared effects related to the FSI of the two final pions. Moreover, in the 
isoscalar case, by comparison with the scalar form factor, we see explicitly 
that this is a universal effect related to the quantum numbers of the 7T7T state. 

The particular dynamics leading to this final state gives rise also to local 
contributions, which are different in each case. We saw in section |3| that 
at the usual reference scale \x = M p the contribution from the local term 
[Lg + 2Ll](M p ) is small. For the weak amplitudes this needs to be fur- 
ther investigated; the usual factorization models ||4f| amount to C|(M p ) = 
Cf(M p )=Cf{M p ) = 0. 

The main effects of the short-distance dynamics, not related to FSI, are 
contained in the particular values of the different amplitudes ai(s) at s = 0. 
This physics needs to be analyzed independently, because it cancels out from 
the Omnes relation. The Omnes factor only allows us to relate the amplitudes 
at two different values of s, but does not give any information on their global 
normalization. 

Taking a low subtraction point where higher-order corrections are ex- 
pected to be small, we can just multiply the tree-level formulae (|4.4j) with 
the experimentally determined Omnes exponentials, as done in ref. |TJ. For 
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/ = we already have the result obtained in the previous section, 



^q(Mk, 0) = 1.55 ± 0.10 , (4.13) 



which improves the lowest-order estimates made in refs. []I],^J] at sq = M 2 . 

In the 1 = 2 channel the inelasticity effect is absent at least up to 1.6 
GeV. Evaluating the once-subtracted dispersive integral over the measured 
phase-shifts up to z = (1.6 GeV) 2 , we get 

^ 1} (M|,0) = 0.92 ±0.03, (4.14) 

to be compared with the earlier estimates ^^(Mj^, M 2 ) = 0.96 (z = 1 GeV 2 ) 
@ and ^\M 2 K ,Ml) = 0.89 ± 0.03 @. The error bar in (|7L|) takes into 



account uncertainties in the fits to the phase-shift data and higher energy 
contributions. 

The corrections induced by FSI in the moduli of the decay amplitudes 
Ai generate an additional enhancement of the A J = 1/2 to AI = 3/2 ratio, 

%{Ml, 0)/3? 2 (M|, 0) = 1.68 ± 0.12 . (4.15) 

This factor multiplies the enhancement already found at short distances. This 
is a quite large correction, which improves previous calculations of Ai(Mg). 
Taking the 3?/ correction into account, the experimental Aj amplitudes imply 
the following corrected values for the lowest-order AS = 1 chiral couplings: 



9&+ g 927 



3.3 , |^27| « 0.31. (4.16) 



These "experimental" numbers are not very far from the short-distance es- 
timates obtained in the first of refs. 



5 Standard Model prediction of e' je 

One further subtlety has to be taken into account in the discussion of CP- 
violating isospin amplitudes. The derivation of eq. fl2.2|) for the absorptive 
parts makes use of Time-Reversal invariance, so that the procedure can be 
strictly applied only to CP-conserving amplitudes. This is not a problem, 
however, because we are working to first order in the weak Fermi coupling. 

The CP-odd phase is hidden in the Wilson coefficients of the short- 
distance AS = 1 Lagrangian (^.1|), which can be decomposed as 

C l {v) = z l {v) + ry l {v) ; r = -^SL . ( 5 .1) 

ud us 
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Since CP violation is only originated by the short-distance ratio of CKM 
matrix elements r, we can always write 



Aj = Af p + TA^ (5.2) 

and apply the Omnes procedure to the amplitudes Af p and Aj , which 
respect Time-Reversal invariance. 

The CP-conserving piece of r Aj is negligible in comparison with Af p . 

Therefore, in a more standard notation, KeAi pa Af p and ImAi = Im(r) A T , 
where "real" and "imaginary" refer to CP-even and CP-odd since the ab- 
sorptive phases have been already factored out through Ai = Aj e lS o . 

The most striking consequence of the correction factors 3?o,2 is a sizeable 
modification of the numerical short-distance estimates^ for the direct CP- 
violation parameter e'/e. A handy way of writing this quantity, used in all 



theoretical short-distance calculations up to date, can be as follows [fL5 



£ 7 =lm(V*V td ) e» [P<va> _ pGW] , (5.3) 

where the phase $ = $ £ / — $ e ~ and the quantities P^ 1 / 2 ) and p( 3 / 2 ) contain 
the contributions from the hadronic matrix elements of four-quark operators 
with AI = 1/2 and 3/2 respectively: 

P a/2) = r <Q»)o (1 - JB ) , 

i 

p(3/2) = - J2Vi( U ) (Qi( U ))2- (5.4) 



Here, (Qi)j = ((inr)i\Qi\K), r and to are given by 



2\e\ ReA ' ReA ' v ' ' 

and the parameter 

1 (lmA 2 )iB f . 

co> lnm 

parametrizes isospin breaking corrections. 

A detailed analysis of e'/e, within the Standard Model, will be given in 
ref. |[20|| . Here we just want to illustrate the important role of FSI and how 



their proper inclusion modifies the e'/ e prediction in a very important way. 

Since the hadronic matrix elements are quite uncertain theoretically, the 
CP-conserving amplitudes KeAj, and thus the factors r and u, are set to 



5 The correction factors ft^ (M% , M%) were already considered in ref. |5^| to estimate 
e'/e within the SU(2) L ® SU(2) R ® U(l) model of CP violation. 
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their experimentally determined values; this automatically includes the FSI 
effect. All the rest in the numerator is theoretically predicted via short- 
distance calculations, because the leading contributions come |TB[ from the 
operators Q 6 and Q 8 whose matrix elements cannot be directly measured 
from K — > 2tt decay rates. 

As a consequence, since the relevant matrix elements (Qe^i are usually 
taken from large-iV^ estimates [15] or lattice calculations which do not 
include FSI corrections, this procedure produces a mismatch with the FSI 
included phenomenologically in the values of r and u. This can be easily 
corrected, introducing in the numerator the dispersion factors 3?j that we 
have estimated. This implies [l|] a large enhancement of the predicted value 
of e'/e by roughly a factor of 2. 

To a very good approximation, the Standard Model prediction for e'/e 
can be written, up to global factors, as [ITH] 



IB 



QABi 3/2) 



(5.7) 



where B^ 2 ^ and B^ 1 *' parametrize the matrix elements of the QCD penguin 
operator Q G and the electroweak penguin operator Q 8 , respectively, in units 
of their vacuum insertion approximation values. These parameters are usu- 
ally taken to be (from 1/Nq considerations |T5j and Lattice calculations |TH] ) 

0.8 ± 0.2, while the isosp in-breaking factor 

With those in- 



,(3/2) 



B 



(1/2) 



is set to Qib 



1.0 ± 0.3 and B, 
0.25 



(3/2) 



56 



57,58 



with large uncertainties 
puts, there is a rather large numerical cancellation between the two terms in 
eq. (|5.7|) , which results in a predicted central value |I|,[l(| e'/e ~ 7.0 x 10~ 4 . 

Since those estimates do not include FSI effects, their values should be 
multiplied by the appropriate factors 3ft/. Notice, that the Omnes procedure 
can be also applied to the individual matrix elements (Qi)j. In order to 
avoid any possible double counting, we will take as the starting point of our 



analysis the large-iVc estimate for the relevant matrix elements |38 

= 1 ; 



B, 



(1/2) 



N C - 



B. 



(3/2) 



1.0 



(5.8) 



FSI only appear at next-to- leading order in the 1/Nc expansion and, there- 
fore, correct the leading values ( |5.8|) . 

The corrected e' / e prediction can be easily obtained, taking into account 
the following points: 

1. The penguin operator Q§ transforms as (8l, 1r) under chiral transfor- 
mations. At lowest order in the chiral expansion, it corresponds to 
the first operator in eq. ( |4.3|) (from eq. ( |5.8j ) one actually gets the Q§ 
contribution to the chiral coupling g 8 , in the large- Nq limit). The FSI 
corrections induced by pion chiral loops modify B^ 2 ^ as follows 



B 



(1/2) 



B, 



(1/2) 



N C - 



Ko(M; 



1,0) 



1.55 
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2. The electroweak penguin operator Q$ corresponds to the chiral operator 
proportional to qew m eq. Q4.3|) . As shown in (|4.4|) , it contributes 



to the two isospin amplitudes, although we only need here the / = 
2 piece (Qs)2- The K — > 2n matrix element is not proportional to 
M\ — M% because the needed SU(3) breaking is provided by the quark 
charge matrix (the chiral operator is identically zero for Q = I). The 
presence or not of this factor does not change the Omnes summation 
(the corresponding zero just factors out whenever is present). One gets 
then 

x sft 2 (M£,0) = 0.92 . 



R (3/2) _ R (3/2) 



N C - 



3. The isosp in-breaking correction coming from ir°-r] mixing has been 
recently calculated at 0(p 4 ) in the chiral expansion, with the result 
Qib — 0.16 ± 0.03 ||59|| . This value is smaller than the previous lowest- 



order estimate VL IB w 0.25 The term B^ 2 ^Q IB in eq. (|5.7|) should 



be multiplied by 9?2 and not by 3?o, because it corresponds to two final 
pions with 1 = 2. Thus, 



R (V2) _ R (l/2) 



n IB x %(a4,o) = o.i5. 



The large FSI correction to the 1 = amplitude gets reinforced by the 
mild suppression of the 1 = 2 contributions. The net effect is a large en- 
hancement of eq. ( [5.71 ), by a factor 2.4, pushing the predicted central value 
from 7.0 x 10~ 4 ||,|l| to 

e7£=17xl0- 4 , (5.9) 

which compares well with the present experimental world average [^l|] in 
eq- O)- 

A more careful analysis, taking into account all hadronic and quark- 
mixing inputs 0] gives the result quoted in eq. ( |1.2| ) for the Standard Model 
prediction of e'/e. 



6 Discussion 



Many attempts have been made to compute the isospin amplitudes Aj from 
first principles ||^, [T5| ; |T6| , p~8| . p~9] . ^Q| . |5^ ? |60] |66[| . Although those calculations have 



provided encouraging results, we are still far from getting accurate predic- 
tions. Nevertheless, a qualitative understanding of the K — > nn transition 
amplitudes is now emerging. 

The strong rescattering of the two final pions generates important cor- 
rections to the kaon decay amplitudes, enhancing the 1 = piece by about 
50% and originating a mild suppression of the 1 = 2 one. FSI alone cannot 
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explain the measured ratio of A J = 1/2 to AI = 3/2 transition amplitudes, 
but they constitute a very important ingredient which reinforces the enhance- 
ment already found at short distances. Combined with the 1/Nc and %PT 
expansions, the calculation of the Omnes factors £Iq 7 2(Mk,0) allows for a 
reliable estimate of e'/e ll|,^0 . 



The lowest-order approximation in the 1/Nc expansion does not provide 
a good starting point to analyze the CP-conserving K — > 2n amplitudes, 
because the anomalous dimensions of the most important operators Qi are 



zero at this order [p4 ] . Thus, at lowest order in 1/Nc one misses the dominant 
physics leading to the well-known short-distance enhancement. That makes 
difficult to perform precise predictions for the K — > 2tt decay rates. 

The situation is different for the CP-violating amplitudes, which are com- 
pletely dominated by Q 6 and Q 8 . These are precisely the only operators 
which have a non-zero anomalous dimension at leading order in the 1/Nc 
expansion. The large-iVc- approximation works rather well for those opera- 
tors |]54|,|60| and their matrix elements can be safely estimated within a 30% 
accuracy, once the large infrared logarithms related to FSI are properly taken 
into account. 

In the large-iVc limit the four-quark operators factorize into currents 
which have a known chiral realization at very low energies. The factorization 
of the operators Qq and Q$ leads to scalar (pseudo-scalar) currents which 
are not directly measurable; their matrix elements are determined with x?T 
techniques at leading (next-to-leading for Q$) non-trivial order in the mo- 
mentum expansion. This fixes the Qq and Q$ contribution to the AS = 1 
%PT couplings in the large-iVc limit fPH, ^0f] ■ A reliable determination of the 
corresponding K — > 2ir transition amplitudes can then be performed at low 
s values where chiral loop corrections are smaller. Once this is accomplished, 
the Omnes dispersive factors allow us to evolve this result to the physical 
s = M\ point, resumming the large chiral corrections associated with FSI. 

The usual vacuum insertion estimate of (Qe)o, adopted in some Standard 
Model calculations of e'/e Hl5| , [li6| , corresponds to the lowest non-trivial order 
in both the 1/Nc and %PT expansions. This naive estimate misses the 



large enhancement generated by one-loop %PT corrections [12-14], which 



originates mainly in the strong rescattering of the two final pions with / = 
J = 0. The FSI correction destroys the accidental numerical cancellation 
between the Q G and Q 8 contributions in eq. ( |5.7|) , producing a large increase 
in the resulting prediction of e'/e. The size of the FSI effect can be already 
determined with the one-loop %PT calculation. The Omnes resummation is 
only needed to perform a reliable estimate of higher-order corrections and 
pin down their associated uncertainties. 

More work is still needed in order to get a precise quantitative description 
of kaon decays. In the meanwhile, our analysis demonstrates that it is at least 
possible to pin down the value of e'/e with an accuracy of about 30%. Within 
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the present uncertainties, the resulting Standard Model theoretical prediction 
is in good agreement with the measured experimental value, without any need 
to invocate a new physics source of CP violation. 
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A Experimental n—n phase shifts 



For the experimental phase shifts of tt-tt scattering with J = and total 
isospin / = or 2 we used a simple parametrization by A. Schenk |67J that 
works in the elastic region. A more involved coupled-channel analysis is 
needed above the first inelastic threshold ||68|| . For generic I and J, the 



parametrization provided in [37] is given by: 



tan Sj(s) 



x 




(A.l) 



The threshold expansion of the scattering amplitude is reproduced by setting 



AM 2 



'J 



(A.2) 



For given J, J there are four parameters: 



a J j, b*j, c J j 



and s 1 ,. The numerical 



values of aj, bj (which parametrize the threshold behaviour of the scattering 
amplitude) have been determined by means of x?T, while the remaining two 
parameters have been extracted from the experimental data. The details of 

Here we only compile the 
= and / = 2, in S-wave, 



the analysis can be found in the original work |67 



numerical values of the parameters. For isospin I 
the values of the threshold parameters are: 



0.20 



b° 



0.24 



-0.042 



hi 



-0.075. 



(A.3) 
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For the other parameters Cj and Sj we have taken the range of values deter- 
mined in |>7 in order to take into account the experimental uncertainties. 



For I 
For I 



0, we used c{] 
2, we used C 2 , 



= 0.008, 0.0, -0.015 and = 840, 865, 890 MeV. 
and s 2 , = -920 2 , -685 2 , -555 2 MeV 2 . The three 



values for each single parameter correspond to the three solid curves shown in 
Figure [I] for each isospin. The central line corresponds to the best fit in 



while the other two extremes enclose the region covered by the experimental 
data considered in p7|. 




200 



400 



600 



800 1000 
E [MeV] 



Figure 1: Phase shifts £o' 2 ( s ) used in our numerical analyses. The solid lines 
enclose the range covered by the experimental data, while the dashed lines 
show the unitarized lowest-order %PT prediction. 

The lowest-order xPT prediction of the phase shifts is also shown with a 
dashed line in Figure [l]. This corresponds to the expression 

tan5 0;2 ( S ) = ^— a^s) (2s - M 2 ; 2M 2 - s) ; ow(«) = ^/l-^, 

(A.4) 

which is a unitarization of the usual prediction, valid at low values 

of s where tan 5 ~ 5. The lowest-order %PT prediction fails already at 
relatively low energies ~ 500 MeV, specially for 1 = 0. In the 1 = case it 
underestimates the experimental phase shift, while in the 1 = 2 case it gives 
a too large, in absolute value, phase shift. 
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B K — ► nn matrix elements at one loop 

The one-loop contribution to the physical CP-conserving K — > tttt isospin 
amplitudes has been computed in Refs. [|l|,[l3|]. In this appendix the same 
amplitudes are calculated at a generic value of the squared invariant mass 
s = (p^i +p n 2) 2 - The complete next-to-leading correction is of order p 4 in 
the chiral expansion and includes one-loop contributions generated by the 
lowest-order p 2 Lagrangian (4.3) and tree-level contributions coming from 



order p 4 counterterms 0,0, 



In the analysis of the CP-conserving amplitudes we have neglected the 
tiny electroweak corrections which are proportional to e 2 gEw & t leading order 
in the chiral expansion. We then decompose the isospin amplitudes as follows 

Ms) = 5i(«) (s - Ml) + 5~ ai {s) (Ml - Ml) , (B.l) 
where Sai(s) is zero at lowest order. In addition, we define clq = + 

~ ~ (%) ~ (27) 

and 5a = 5%' + 5% , thus explicitly separating the octet and 27-plet 
contributions to the 1 = amplitude. 

At 0(p 4 ) the octet 1 = function takes the form 



~(8) 



^V ud V: s V2Ug 8 jl - \ (25 - Ml) B(M*,M 2 ,s) 
+ Us-M 2 K + Ml) B(M 2 K , Ml s) + \Ml B(M 2 , M 2 , s) 



+ 



+ 



{Ml - Ml) 

12M 2 
67 1 
12^+6 
1 



6 

(3M 2 K - AMI) B(M 2 K , Ml, Ml) + M 2 K B(M 2 K , M 2 , M 2 n 
7 {Ml- Ml) , M 2 



2M 2 K 



P 



36 r "'' ' 2M 2 
-16(2M| + Ml) Ll(fi) - A(M 2 K + 3M 2 n ) L^) + M 2 K C*(n) 
M 4 

+ mIc 8 m- j ^c 8 m + sC 8 m 



(B.2) 



while 5a^ is given by 
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+ 2 m ) - _ ^ + - j . (B.3) 



M 2 K -M 2 \2 18^ ^7 2M 2 K (M 2 K - M 2 ) 
At the one-loop level, the 27-plet 7 = function is 



4 27) = 



'V2 



V ud V: s V2fJ- g 27 \l - \ (2s - M 2 ) B(M 2 , M 2 , s) 



+ - (s - M 2 K + M 2 ) B(M 2 K , M 2 K , s) - -Ml B(M 2 , M 2 , s) 



+ 



{Ml - Ml) 



12M 2 

TV 



(3M| - AMI) B(Ml, Ml, Ml) - AM 2 K B(M 2 K , M 2 , Ml) 



+ 



67 

il 
i 
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+ 12 ^ n — 6~ ^ K ~ A ^ v + 



(M: 



1-M 



2M2 



— (/^ + Afi v - 5n K ) 



2 ^/i^ //„-J 



2M£ 



/ 2 



-16(2M 2 + M 2 ) LJ(ii) - 4(M^ + 3M 2 ) L r M + M'l Cf^) 



+MlCf{») 



M: 



Ml 



* C 27 (pi) + sC?(ri 



(B-4) 



and 



5a, 



^(27) 



G F 



^V ud V: s V2f n -g 27 I -Ml [AB(M 2 ,M 2 ,s)-3B(M 2 K ,M 2 K ,s) 



- ^(5M| + AMI) B(M 2 K , Ml, M 2 J - l -M 2 K B(M 2 K , M 2 , M 2 n ) 



+ 



P 



M, 



Ml Cf(/x) + M 2 K CfQi) + Ci 7 (n) 



Ml 



M 2 



Ml - Ml 



— + -fi v - 2/i^ ) + 



Mi 



5 n 13 
- g ^ - 2 + — /i X 



(/^ - /^) \ ■ (B.5) 



Finally, the I = 2 a 2 function is 



a 2 = ~^V ud V: s y/ w ^ |l + I ( a - 2M 2 ) 5(M 2 , M 2 , S ) 
(M£-M 2 



+ 



24M2 

7T 



23 
15 
1 



(15M| - 8M 2 ) S(M|, M 2 , Ml) + M 2 K B(M 2 K , M 2 , M 2 n ] 
7 



(M 2 K - Ml) . 

^2 ^ ~ 4 - g + 4^ ( 5 /^ - - 4//k 



16(2M| + M 2 )LM - 4(M 2 + SM 2 )^) + M* C?{n) 
+ MlCf^) + sCf{j), 



(B.6) 
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while <5a 2 is given by 



+ 



i 



*a> = -^W s {-^ [M^f (,x) + M|Cf(/i) 
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1 



M£ £(M£, M 2 , M 2 ) - (M 2 K + 8M^) B(M^, M 2 , M*) 

Ml ^ 



-(^ + 3^ + 4^) + Mi _ M , 



(B.7) 



The one- loop function B(Mf,Mf ,p ) is defined as follows 



B(Ml,Miy 

+ ^ I 



16vr 2 / 2 
M 2 M| 

p2 ' ^2 



-1 + In -f 



In 



- Ml 



Mf + A 1/2 (p 2 ,M?,M%) 



f - M\ - Ml - A 1 / 2 (p2 a M 2 , Mf) 
for p 2 > (Mi + M 2 ) 2 and p 2 < (Mj - M 2 ) 2 , while 



(B.8) 



5(M 2 ,M 2 2 ,^ 



lQn 2 p 



-1 + ln 



'Mf 



1 /Mf 



1 + 



M 2 M| 

p2 p2 



-A 1 



Ml Mf\ 
p 2 , p 2 I 



arctan 



(p*-Ml-Mi)J-\(p 2 ,Ml,M 2 ) 



(p 2 - Ml - Ml) 2 - 2MIMI 



(B.9) 



for (M x - M 2 ) 2 < p 2 < (Mi + M 2 ) 2 and -tt/2 < arctan (x) < tt/2. The 
function A(x, y, z) is given by 



X(x,y,z) = (x + y — z) 2 — Axy . 



(B.10) 



The function _B(M 2 ,Mf p 2 ) is related to the function J{s) introduced in 
ref. in the following way 



f 2 B(Mf 1 M 2 1 p 2 ) = -J(s) 



For Mi = M 2 = m one gets 



1 / M 2 M 2 , M?\ 

In-f + TT^T^ln-^ . (B.ll) 



16vr 2 I /i 2 M 2 - Mf M| 



M 2 



/ 2 B(M 2 ,MV) = -J(s) + — (h^ + I 
where J(s) for Mi = M 2 has been given in eq. ( p.5|) . 



(B.12) 
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The parameters \xp (P = it, K, 77) contain a logarithmic dependence on 
the chiral renormalization scale \l generated by one-loop corrections. They 
are defined as 

M 2 M 2 

< B13 > 

The explicit /i-dependence of \ip and the functions B(Mf, M 2 , p 2 ) is canceled 
by the local contributions L\ (i = 4, 5), Cf (i = 1, . . . 6), Cf (z = 1, . . . 6) 
andCf(i = l,...5). 



C Comments on recent literature 

The ideas put forward in our first letter have been further discussed in 
several recent papers by other authors. We would like to make here some 
brief comments on these works. 



C.l 



34j| that the Omnes exponential depends 
that reference the FSI enhancement 



It has been pointed out in ref. 
on the chosen subtraction point. In 
is minimized, by taking the highest subtraction point below the physical 
cut, s = 4:M 2 . This trivial fact is then used to argue that our evaluation 
of FSI corrections is unreliable. Our detailed analysis of the subtraction 
point dependence in section [5] shows that this claim is unfounded. Taking 
a higher value of Sq one is just shifting FSI corrections from the Omnes 
exponential to the amplitude in front, but the physical result is of course the 
same. At s — 4M^ there is a large one-loop correction to the amplitude, 
which has been overlooked in ref. EM. Moreover, s = 4M^ is a bad choice 



of subtraction point, because the corresponding Taylor expansion has zero 



convergence radius [35 



Ref. states that it is not precisely known at which value of s existing 
lattice estimates correspond to. It is suggested that future lattice calculations 
could obtain the K — > nir weak matrix elements at threshold (sq = 4M 2 ) 
and the 9?j(Mg-, s ) correction factors could then be used to get the physical 
amplitudes. 

The s ambiguity mentioned in ref. 34] is not present in the low-energy 
chiral expansion. Our e'/e calculation is based on a large- Nq evaluation of 
the couplings of the AS* = 1 Lagrangian. Once these chiral couplings 
are determined the K — > amplitudes can in principle be computed at 
any value of s. Higher-order chiral corrections are of course smaller at lower 
values of momenta, which makes advisable to use the chiral expansion al low 
s values. 
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C.2 

At lowest order in xPT, the four-quark operator Q 8 induces the 0(p°) chiral 
term proportional to the coupling Qew- The corresponding one-loop cor- 
rection has been analyzed in ref. pE5|] , where a small positive contribution is 



obtained. This result agrees with a recent dispersive calculation of B^ 2 ^ [BS 



which finds Bf /2) = 1.11 ± 0.16 ± 0.23. 

FSI generate instead a small suppression of the 1 = 2 amplitude. How- 
ever, there are other chiral corrections not related to FSI which appear at the 
one-loop level; they are included in the value of a 2 (0). Since for / = 2 the FSI 
effect is small, other correction could be equally important and even reverse 
the sign of the correction for the physical amplitude G^Mj^). A detailed 
one-loop analysis will be presented elsewhere pT7|| . 

C.3 

A simplified (the dispersive integral over the phase-shift is not exponentiated) 
version of our Jij- factors has been used in ref. JF(J , which advocates a different 



and conceptually incorrect interpretation of the chiral corrections related to 
FSI. 

In this reference a non-subtracted dispersion relation is used. The re- 
sulting divergence in Iq ^s) is regulated cutting the dispersive integral at 
the upper edge z, and making the ad- hoc identification z = u 2 , with v the 
short-distance scale governing the Wilson coefficients of the effective four- 
quark Hamiltonian (|4.1|) . This generates a v dependence in the dispersive 
integral which is claimed to cancel the renormalization-scale dependence of 
the Wilson coefficients. 

The choice z = v 2 and the associated identification of infrared and ul- 
traviolet logarithms is arbitrary and cannot be correct. The FSI logarithms 
have nothing to do with the underlying short-distance physics. The Omnes 
factor Qj(s,s ) relates the isospin amplitudes at two different points s and 
s, but is unable to fix the global normalization. The short-distance informa- 
tion is hidden in ar(0) which, moreover, is independent of the scale v. Thus, 
the cancellation of v dependences must be accomplished even in the absence 
of FSI. 

The argument can be better seen analyzing the scalar form factor, which 
has the same FSI logarithms but a different short-distance contribution. In 
fact, instead of working with hadronic matrix elements of the scalar current, 
we can take the corresponding matrix elements of the divergence of the asso- 
ciated vector current. Both quantities are trivially related by a quark-mass 
factor, through a Ward identity. The FSI phenomena and the associated %PT 
logarithms are of course identical; but now, there is no anomalous dimension. 
Since there is no short-distance renormalization scale z/, the identification 
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z = v 2 is then meaningless. 

The same argumentation can be repeated with the Omnes summation of 
FSI effects in the pion vector form factor |25| . Again, this is a renormalizat ion- 
group invariant quantity (no anomalous dimension) and does not make any 
sense to identify FSI phenomena with non-existing short-distance logarithms. 
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